Genome-wide CRISPR screens identify the YAP/TEAD axis as a driver of persister cells in EGFR mutant lung cancer.

Most lung cancer patients with metastatic cancer eventually relapse with drug-resistant disease following treatment and EGFR mutant lung cancer is no exception. Genome-wide CRISPR screens, to either knock out or overexpress all protein-coding genes in cancer cell lines, revealed the landscape of pathways that cause resistance to the EGFR inhibitors osimertinib or gefitinib in EGFR mutant lung cancer. Among the most recurrent resistance genes were those that regulate the Hippo pathway. Following osimertinib treatment a subpopulation of cancer cells are able to survive and over time develop stable resistance. These 'persister' cells can exploit non-genetic (transcriptional) programs that enable cancer cells to survive drug treatment. Using genetic and pharmacologic tools we identified Hippo signalling as an important non-genetic mechanism of cell survival following osimertinib treatment. Further, we show that combinatorial targeting of the Hippo pathway and EGFR is highly effective in EGFR mutant lung cancer cells and patient-derived organoids, suggesting a new therapeutic strategy for EGFR mutant lung cancer patients.

Overcoming osimertinib resistance with AKT inhibition in EGFRm-driven Non-Small-Cell-Lung-Cancer with PIK3CA/PTEN alterations.

PURPOSE: Osimertinib is an epidermal growth factor receptor tyrosine kinase inhibitor (EGFR-TKI) indicated for the treatment of EGFR mutated (EGFRm)-driven lung adenocarcinomas. Osimertinib significantly improves progression-free survival in first-line treated patients with EGFRm advanced NSCLC. Despite the durable disease control, the majority of patients receiving osimertinib eventually develop disease progression. EXPERIMENTAL DESIGN: ctDNA profiling analysis on-progression plasma samples from patients treated with osimertinib in both first (Phase 3, FLAURA trial) and second-line trials (Phase 3, AURA3 trial) revealed a high prevalence of PIK3CA/AKT/PTEN alterations. In vitro and in vivo evidence using CRISPR engineered NSCLC cell lines and PXD models support a functional role for PIK3CA and PTEN mutations in the development of osimertinib resistance. RESULTS: These alterations are functionally relevant as EGFRm NSCLC cells with engineered PIK3CA/AKT/PTEN alterations develop resistance to osimertinib and can be re-sensitized by treatment with the combination of osimertinib and the AKT inhibitor capivasertib. Moreover, xenograft and PDX in vivo models with PIK3CA/AKT/PTEN alterations display limited sensitivity to osimertinib relative to models without alteration, and in these double mutant models capivasertib and osimertinib combination elicits an improved anti-tumor effect versus osimertinib alone. CONCLUSIONS: Together, this approach offers a potential treatment strategy for patients with EGFRm-driven NSCLC that have a sub-optimal response, or develop resistance, to osimertinib through PIK3CA/AKT/PTEN alterations.

Modeling Flexible Protein Structure With AlphaFold2 and Crosslinking Mass Spectrometry.

We propose a pipeline that combines AlphaFold2 (AF2) and crosslinking mass spectrometry (XL-MS) to model the structure of proteins with multiple conformations. The pipeline consists of two main steps: ensemble generation using AF2 and conformer selection using XL-MS data. For conformer selection, we developed two scores-the monolink probability score (MP) and the crosslink probability score (XLP)-both of which are based on residue depth from the protein surface. We benchmarked MP and XLP on a large dataset of decoy protein structures and showed that our scores outperform previously developed scores. We then tested our methodology on three proteins having an open and closed conformation in the Protein Data Bank: Complement component 3 (C3), luciferase, and glutamine-binding periplasmic protein, first generating ensembles using AF2, which were then screened for the open and closed conformations using experimental XL-MS data. In five out of six cases, the most accurate model within the AF2 ensembles-or a conformation within 1 Å of this model-was identified using crosslinks, as assessed through the XLP score. In the remaining case, only the monolinks (assessed through the MP score) successfully identified the open conformation of glutamine-binding periplasmic protein, and these results were further improved by including the "occupancy" of the monolinks. This serves as a compelling proof-of-concept for the effectiveness of monolinks. In contrast, the AF2 assessment score was only able to identify the most accurate conformation in two out of six cases. Our results highlight the complementarity of AF2 with experimental methods like XL-MS, with the MP and XLP scores providing reliable metrics to assess the quality of the predicted models. The MP and XLP scoring functions mentioned above are available at https://gitlab.com/topf-lab/xlms-tools.

Knowledge graph-based recommendation framework identifies drivers of resistance in EGFR mutant non-small cell lung cancer.

Resistance to EGFR inhibitors (EGFRi) presents a major obstacle in treating non-small cell lung cancer (NSCLC). One of the most exciting new ways to find potential resistance markers involves running functional genetic screens, such as CRISPR, followed by manual triage of significantly enriched genes. This triage process to identify 'high value' hits resulting from the CRISPR screen involves manual curation that requires specialized knowledge and can take even experts several months to comprehensively complete. To find key drivers of resistance faster we build a recommendation system on top of a heterogeneous biomedical knowledge graph integrating pre-clinical, clinical, and literature evidence. The recommender system ranks genes based on trade-offs between diverse types of evidence linking them to potential mechanisms of EGFRi resistance. This unbiased approach identifies 57 resistance markers from >3,000 genes, reducing hit identification time from months to minutes. In addition to reproducing known resistance markers, our method identifies previously unexplored resistance mechanisms that we prospectively validate.

Pharmaceutical Reactivation of Attenuated Apoptotic Pathways Leads to Elimination of Osimertinib Drug-Tolerant Cells.

Osimertinib is an EGFR tyrosine kinase inhibitor (TKI) with proven clinical efficacy; however, acquired resistance presents an obstacle to curing EGFR-driven disease. Recent studies have shown that drug-tolerant persister cells (DTP) have a distinct transcriptional profile that may confer specific vulnerabilities. By definition these cells avoid apoptosis, yet little is known about how their survival is regulated. We found that paradoxically, the proapoptotic gene BIM was upregulated in osimertinib DTPs, and cotreatment with BH3 mimetics could trigger DTP cell death. Furthermore, cIAP proteins, antiapoptotic members of the extrinsic pathway, were significantly elevated in DTPs. cIAP antagonists could block DTP formation as an up-front combination, and could eliminate preformed DTPs. Critically, when treated at the time of maximal osimertinib response, cIAP or MCL1 inhibitor treatment could significantly attenuate the regrowth of EGFRm cell line mouse xenografts. Finally, we show that apoptosis can be maximized in cell lines with acquired osimertinib resistance by combining BH3 or SMAC mimetics with agents that target the resistance driver in these models. Taken together, these data suggest novel therapeutic strategies at the point of minimal residual disease or full osimertinib resistance for patients in this critical area of unmet need. Significance: These studies uncover strategies to use targeted agents that activate apoptosis in non-small cell lung cancer cells that survive initial EGFR TKI treatment.

A YAP/FOXM1 axis mediates EMT-associated EGFR inhibitor resistance and increased expression of spindle assembly checkpoint components.

Acquired resistance to tyrosine kinase inhibitors (TKIs) of epidermal growth factor receptor (EGFR) remains a clinical challenge. Especially challenging are cases in which resistance emerges through EGFR-independent mechanisms, such as through pathways that promote epithelial-to-mesenchymal transition (EMT). Through an integrated transcriptomic, proteomic, and drug screening approach, we identified activation of the yes-associated protein (YAP) and forkhead box protein M1 (FOXM1) axis as a driver of EMT-associated EGFR TKI resistance. EGFR inhibitor resistance was associated with broad multidrug resistance that extended across multiple chemotherapeutic and targeted agents, consistent with the difficulty of effectively treating resistant disease. EGFR TKI-resistant cells displayed increased abundance of spindle assembly checkpoint (SAC) proteins, including polo-like kinase 1 (PLK1), Aurora kinases, survivin, and kinesin spindle protein (KSP). Moreover, EGFR TKI-resistant cells exhibited vulnerability to SAC inhibitors. Increased activation of the YAP/FOXM1 axis mediated an increase in the abundance of SAC components in resistant cells. The clinical relevance of these finding was indicated by evaluation of specimens from patients with EGFR mutant lung cancer, which showed that high FOXM1 expression correlated with expression of genes encoding SAC proteins and was associated with a worse clinical outcome. These data revealed the YAP/FOXM1 axis as a central regulator of EMT-associated EGFR TKI resistance and that this pathway, along with SAC components, are therapeutic vulnerabilities for targeting this multidrug-resistant phenotype.

SSB1/SSB2 Proteins Safeguard B Cell Development by Protecting the Genomes of B Cell Precursors.

Induction of programmed DNA damage and its recognition and repair are fundamental for B cell development. The ssDNA-binding protein SSB1 has been described in human cells as essential for the recognition and repair of DNA damage. To study its relevance for B cells, we recently developed Ssb1 -/- and conditional Ssb1 -/- mice. Although SSB1 loss did not affect B cell development, Ssb1 -/- cells exhibited compensatory expression of its homolog SSB2. We have now generated Ssb2 -/- mice and show in this study that SSB2 is also dispensable for B cell development and DNA damage response activation. In contrast to the single loss of Ssb1 or Ssb2, however, combined SSB1/2 deficiency caused a defect in early B cell development. We relate this to the sensitivity of B cell precursors as mature B cells largely tolerated their loss. Toxicity of combined genetic SSB1/2 loss can be rescued by ectopic expression of either SSB1 or SSB2, mimicked by expression of SSB1 ssDNA-binding mutants, and attenuated by BCL2-mediated suppression of apoptosis. SSB1/2 loss in B cell precursors further caused increased exposure of ssDNA associated with disruption of genome fragile sites, inefficient cell cycle progression, and increased DNA damage if apoptosis is suppressed. As such, our results establish SSB1/2 as safeguards of B cell development and unveil their differential requirement in immature and mature B lymphocytes.

The Systems Architecture of Molecular Memory in Poplar after Abiotic Stress.

Throughout the temperate zones, plants face combined drought and heat spells in increasing frequency and intensity. Here, we compared periodic (intermittent, i.e., high-frequency) versus chronic (continuous, i.e., high-intensity) drought-heat stress scenarios in gray poplar (Populus× canescens) plants for phenotypic and transcriptomic effects during stress and after recovery. Photosynthetic productivity after stress recovery exceeded the performance of poplar trees without stress experience. We analyzed the molecular basis of this stress-related memory phenotype and investigated gene expression responses across five major tree compartments including organs and wood tissues. For each of these tissue samples, transcriptomic changes induced by the two stress scenarios were highly similar during the stress phase but strikingly divergent after recovery. Characteristic molecular response patterns were found across tissues but involved different genes in each tissue. Only a small fraction of genes showed similar stress and recovery expression profiles across all tissues, including type 2C protein phosphatases, the LATE EMBRYOGENESIS ABUNDANT PROTEIN4-5 genes, and homologs of the Arabidopsis (Arabidopsis thaliana) transcription factor HOMEOBOX7. Analysis of the predicted transcription factor regulatory networks for these genes suggested that a complex interplay of common and tissue-specific components contributes to the coordination of post-recovery responses to stress in woody plants.

Optical and physical mapping with local finishing enables megabase-scale resolution of agronomically important regions in the wheat genome.

BACKGROUND: Numerous scaffold-level sequences for wheat are now being released and, in this context, we report on a strategy for improving the overall assembly to a level comparable to that of the human genome. RESULTS: Using chromosome 7A of wheat as a model, sequence-finished megabase-scale sections of this chromosome were established by combining a new independent assembly using a bacterial artificial chromosome (BAC)-based physical map, BAC pool paired-end sequencing, chromosome-arm-specific mate-pair sequencing and Bionano optical mapping with the International Wheat Genome Sequencing Consortium RefSeq v1.0 sequence and its underlying raw data. The combined assembly results in 18 super-scaffolds across the chromosome. The value of finished genome regions is demonstrated for two approximately 2.5 Mb regions associated with yield and the grain quality phenotype of fructan carbohydrate grain levels. In addition, the 50 Mb centromere region analysis incorporates cytological data highlighting the importance of non-sequence data in the assembly of this complex genome region. CONCLUSIONS: Sufficient genome sequence information is shown to now be available for the wheat community to produce sequence-finished releases of each chromosome of the reference genome. The high-level completion identified that an array of seven fructosyl transferase genes underpins grain quality and that yield attributes are affected by five F-box-only-protein-ubiquitin ligase domain and four root-specific lipid transfer domain genes. The completed sequence also includes the centromere.

Computational Design of Epitope-Enriched HIV-1 Gag Antigens with Preserved Structure and Function for Induction of Broad CD8+ T Cell Responses.

The partially protective phenotype observed in HIV-infected long-term-non-progressors is often associated with certain HLA alleles, thus indicating that cytotoxic T lymphocyte (CTL) responses play a crucial role in combating virus replication. However, both the vast variability of HIV and the HLA diversity impose a challenge on elicitation of broad and effective CTL responses. Therefore, we conceived an algorithm for the enrichment of CD8+ T cell epitopes in HIV's Gag protein, respecting functional preservation to enable cross-presentation. Experimentally identified epitopes were compared to a Gag reference sequence. Amino-acid-substitutions (AAS) were assessed for their impact on Gag's budding-function using a trained classifier that considers structural models and sequence conservation. Experimental assessment of Gag-variants harboring selected AAS demonstrated an apparent classifier-precision of 100%. Compatible epitopes were assigned an immunological score that incorporates features such as conservation or HLA-association in a user-defined weighted manner. Using a genetic algorithm, the epitopes were incorporated in an iterative manner into novel T-cell-epitope-enriched Gag sequences (TeeGag). Computational evaluation showed that these antigen candidates harbor a higher fraction of epitopes with higher score as compared to natural Gag isolates and other artificial antigen designs. Thus, these designer sequences qualify as next-generation antigen candidates for induction of broader CTL responses.

Sensitivity to PI3K and AKT inhibitors is mediated by divergent molecular mechanisms in subtypes of DLBCL.

Activated B-cell-like (ABC) and germinal center B-cell-like diffuse large B-cell lymphoma (DLBCL) represent the 2 major molecular DLBCL subtypes. They are characterized by differences in clinical course and by divergent addiction to oncogenic pathways. To determine activity of novel compounds in these 2 subtypes, we conducted an unbiased pharmacologic in vitro screen. The phosphatidylinositol-3-kinase (PI3K) α/δ (PI3Kα/δ) inhibitor AZD8835 showed marked potency in ABC DLBCL models, whereas the protein kinase B (AKT) inhibitor AZD5363 induced apoptosis in PTEN-deficient DLBCLs irrespective of their molecular subtype. These in vitro results were confirmed in various cell line xenograft and patient-derived xenograft mouse models in vivo. Treatment with AZD8835 induced inhibition of nuclear factor κB signaling, prompting us to combine AZD8835 with the Bruton's tyrosine kinase inhibitor ibrutinib. This combination was synergistic and effective both in vitro and in vivo. In contrast, the AKT inhibitor AZD5363 was effective in PTEN-deficient DLBCLs through downregulation of the oncogenic transcription factor MYC. Collectively, our data suggest that patients should be stratified according to their oncogenic dependencies when treated with PI3K and AKT inhibitors.

Determination of moderately polar arsenolipids and mercury speciation in freshwater fish of the River Elbe (Saxony, Germany).

Arsenic and mercury are frequent contaminants in the environment and care must be taken to limit their entrance into the food chain. The toxicity of both elements strongly depends upon their speciation. Total amounts of As and Hg as well as their species were analyzed in muscle and liver of 26 fishes of seven freshwater fish species caught in the River Elbe. The median concentrations of As were 162 µg kg(-1) w.w. in liver and 92 µg kg(-1) w.w. in muscle. The median concentrations of total Hg were 241 µg kg(-1) w.w. in liver and 256 µg kg(-1) w.w. in muscle. While this level of Hg contamination of the freshwater fish in the River Elbe is significantly lower than 20 years ago, it exceeds the recommended environmental quality standard of 20 µg Hg kg(-1) w.w. by a factor of 5-50. However, the European maximum level of 500 µg Hg kg(-1) for fish for human consumption is rarely exceeded. Arsenic-containing fatty acids and hydrocarbons were determined and partially identified in methanolic extracts of the fish by HPLC coupled in parallel to ICP-MS (element specific detection) and ESI-Q-TOF-MS (molecular structure detection). While arsenobetaine was the dominant As species in the fish, six arsenolipids were detected and identified in the extracts of liver tissue in common bream (Abramis brama), ide (Leuciscus idus), asp (Aspius aspius) and northern pike (Esox lucius). Four arsenic-containing fatty acids (AsFA) and two arsenic-containing hydrocarbons (AsHC) are reported in freshwater fish for the first time. With respect to mercury the more toxic MeHg(+) was the major species in muscle tissue (>90% of total Hg) while in liver Hg(2+) and MeHg(+) were of equal importance. The results show the high relevance of element speciation in addition to the determination of total element concentrations to correctly assess the burden of these two elements in fish.

PGSB/MIPS Plant Genome Information Resources and Concepts for the Analysis of Complex Grass Genomes.

PGSB (Plant Genome and Systems Biology; formerly MIPS-Munich Institute for Protein Sequences) has been involved in developing, implementing and maintaining plant genome databases for more than a decade. Genome databases and analysis resources have focused on individual genomes and aim to provide flexible and maintainable datasets for model plant genomes as a backbone against which experimental data, e.g., from high-throughput functional genomics, can be organized and analyzed. In addition, genomes from both model and crop plants form a scaffold for comparative genomics, assisted by specialized tools such as the CrowsNest viewer to explore conserved gene order (synteny) between related species on macro- and micro-levels.The genomes of many economically important Triticeae plants such as wheat, barley, and rye present a great challenge for sequence assembly and bioinformatic analysis due to their enormous complexity and large genome size. Novel concepts and strategies have been developed to deal with these difficulties and have been applied to the genomes of wheat, barley, rye, and other cereals. This includes the GenomeZipper concept, reference-guided exome assembly, and "chromosome genomics" based on flow cytometry sorted chromosomes.

Modulation of Ambient Temperature-Dependent Flowering in Arabidopsis thaliana by Natural Variation of FLOWERING LOCUS M.

Plants integrate seasonal cues such as temperature and day length to optimally adjust their flowering time to the environment. Compared to the control of flowering before and after winter by the vernalization and day length pathways, mechanisms that delay or promote flowering during a transient cool or warm period, especially during spring, are less well understood. Due to global warming, understanding this ambient temperature pathway has gained increasing importance. In Arabidopsis thaliana, FLOWERING LOCUS M (FLM) is a critical flowering regulator of the ambient temperature pathway. FLM is alternatively spliced in a temperature-dependent manner and the two predominant splice variants, FLM-ß and FLM-δ, can repress and activate flowering in the genetic background of the A. thaliana reference accession Columbia-0. The relevance of this regulatory mechanism for the environmental adaptation across the entire range of the species is, however, unknown. Here, we identify insertion polymorphisms in the first intron of FLM as causative for accelerated flowering in many natural A. thaliana accessions, especially in cool (15°C) temperatures. We present evidence for a potential adaptive role of this structural variation and link it specifically to changes in the abundance of FLM-ß. Our results may allow predicting flowering in response to ambient temperatures in the Brassicaceae.

Joint Transcriptomic and Metabolomic Analyses Reveal Changes in the Primary Metabolism and Imbalances in the Subgenome Orchestration in the Bread Wheat Molecular Response to Fusarium graminearum.

Fusarium head blight is a prevalent disease of bread wheat (Triticum aestivum L.), which leads to considerable losses in yield and quality. Quantitative resistance to the causative fungus Fusarium graminearum is poorly understood. We integrated transcriptomics and metabolomics data to dissect the molecular response to the fungus and its main virulence factor, the toxin deoxynivalenol in near-isogenic lines segregating for two resistance quantitative trait loci, Fhb1 and Qfhs.ifa-5A. The data sets portrait rearrangements in the primary metabolism and the translational machinery to counter the fungus and the effects of the toxin and highlight distinct changes in the metabolism of glutamate in lines carrying Qfhs.ifa-5A. These observations are possibly due to the activity of two amino acid permeases located in the quantitative trait locus confidence interval, which may contribute to increased pathogen endurance. Mapping to the highly resolved region of Fhb1 reduced the list of candidates to few genes that are specifically expressed in presence of the quantitative trait loci and in response to the pathogen, which include a receptor-like protein kinase, a protein kinase, and an E3 ubiquitin-protein ligase. On a genome-scale level, the individual subgenomes of hexaploid wheat contribute differentially to defense. In particular, the D subgenome exhibited a pronounced response to the pathogen and contributed significantly to the overall defense response.

A synteny-based draft genome sequence of the forage grass Lolium perenne.

Here we report the draft genome sequence of perennial ryegrass (Lolium perenne), an economically important forage and turf grass species that is widely cultivated in temperate regions worldwide. It is classified along with wheat, barley, oats and Brachypodium distachyon in the Pooideae sub-family of the grass family (Poaceae). Transcriptome data was used to identify 28,455 gene models, and we utilized macro-co-linearity between perennial ryegrass and barley, and synteny within the grass family, to establish a synteny-based linear gene order. The gametophytic self-incompatibility mechanism enables the pistil of a plant to reject self-pollen and therefore promote out-crossing. We have used the sequence assembly to characterize transcriptional changes in the stigma during pollination with both compatible and incompatible pollen. Characterization of the pollen transcriptome identified homologs to pollen allergens from a range of species, many of which were expressed to very high levels in mature pollen grains, and are potentially involved in the self-incompatibility mechanism. The genome sequence provides a valuable resource for future breeding efforts based on genomic prediction, and will accelerate the development of new varieties for more productive grasslands.

chromoWIZ: a web tool to query and visualize chromosome-anchored genes from cereal and model genomes.

BACKGROUND: Over the last years reference genome sequences of several economically and scientifically important cereals and model plants became available. Despite the agricultural significance of these crops only a small number of tools exist that allow users to inspect and visualize the genomic position of genes of interest in an interactive manner. DESCRIPTION: We present chromoWIZ, a web tool that allows visualizing the genomic positions of relevant genes and comparing these data between different plant genomes. Genes can be queried using gene identifiers, functional annotations, or sequence homology in four grass species (Triticum aestivum, Hordeum vulgare, Brachypodium distachyon, Oryza sativa). The distribution of the anchored genes is visualized along the chromosomes by using heat maps. Custom gene expression measurements, differential expression information, and gene-to-group mappings can be uploaded and can be used for further filtering. CONCLUSIONS: This tool is mainly designed for breeders and plant researchers, who are interested in the location and the distribution of candidate genes as well as in the syntenic relationships between different grass species. chromoWIZ is freely available and online accessible at http://mips.helmholtz-muenchen.de/plant/chromoWIZ/index.jsp.

Genome interplay in the grain transcriptome of hexaploid bread wheat.

Allohexaploid bread wheat (Triticum aestivum L.) provides approximately 20% of calories consumed by humans. Lack of genome sequence for the three homeologous and highly similar bread wheat genomes (A, B, and D) has impeded expression analysis of the grain transcriptome. We used previously unknown genome information to analyze the cell type-specific expression of homeologous genes in the developing wheat grain and identified distinct co-expression clusters reflecting the spatiotemporal progression during endosperm development. We observed no global but cell type- and stage-dependent genome dominance, organization of the wheat genome into transcriptionally active chromosomal regions, and asymmetric expression in gene families related to baking quality. Our findings give insight into the transcriptional dynamics and genome interplay among individual grain cell types in a polyploid cereal genome.

Ancient hybridizations among the ancestral genomes of bread wheat.

The allohexaploid bread wheat genome consists of three closely related subgenomes (A, B, and D), but a clear understanding of their phylogenetic history has been lacking. We used genome assemblies of bread wheat and five diploid relatives to analyze genome-wide samples of gene trees, as well as to estimate evolutionary relatedness and divergence times. We show that the A and B genomes diverged from a common ancestor ~7 million years ago and that these genomes gave rise to the D genome through homoploid hybrid speciation 1 to 2 million years later. Our findings imply that the present-day bread wheat genome is a product of multiple rounds of hybrid speciation (homoploid and polyploid) and lay the foundation for a new framework for understanding the wheat genome as a multilevel phylogenetic mosaic.

Ragweed (Ambrosia artemisiifolia) pollen allergenicity: SuperSAGE transcriptomic analysis upon elevated CO2 and drought stress.

BACKGROUND: Pollen of common ragweed (Ambrosia artemisiifolia) is a main cause of allergic diseases in Northern America. The weed has recently become spreading as a neophyte in Europe, while climate change may also affect the growth of the plant and additionally may also influence pollen allergenicity. To gain better insight in the molecular mechanisms in the development of ragweed pollen and its allergenic proteins under global change scenarios, we generated SuperSAGE libraries to identify differentially expressed transcripts. RESULTS: Ragweed plants were grown in a greenhouse under 380 ppm CO2 and under elevated level of CO2 (700 ppm). In addition, drought experiments under both CO2 concentrations were performed. The pollen viability was not altered under elevated CO2, whereas drought stress decreased its viability. Increased levels of individual flavonoid metabolites were found under elevated CO2 and/or drought. Total RNA was isolated from ragweed pollen, exposed to the four mentioned scenarios and four SuperSAGE libraries were constructed. The library dataset included 236,942 unique sequences, showing overlapping as well as clear differently expressed sequence tags (ESTs). The analysis targeted ESTs known in Ambrosia, as well as in pollen of other plants. Among the identified ESTs, those encoding allergenic ragweed proteins (Amb a) increased under elevated CO2 and drought stress. In addition, ESTs encoding allergenic proteins in other plants were also identified. CONCLUSIONS: The analysis of changes in the transcriptome of ragweed pollen upon CO2 and drought stress using SuperSAGE indicates that under global change scenarios the pollen transcriptome was altered, and impacts the allergenic potential of ragweed pollen.